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Abstract 

We consider a supersymmetric model where the neutrino mass matrix arises from bilinear and 
trilinear R-parity violation, both restricted by a Z3 flavor symmetry. Assuming flavor blind soft 
supersymmetry (SUSY) breaking conditions, corrected at low energies due to running effects, we 
obtain a neutrino mass matrix in agreement with oscillation data. In particular, a large #13 angle 
can be easily accommodated. We also find a correlation between the reactor and atmospheric 
mixing angles. This leads in some scenarios to a clear deviation from #23 = vr/4. The Lightest 
Supersymmetric Particle (LSP) decay, dominated by the trilinear couplings, provides a direct way 
to test the model at colliders. 

Keywords: supersymmetry; neutrino masses and mixing 



* Electronic address: epeinado@ific.uv.es 



^Electronic address: |avelino.vicente@th.u-psud.fr 



I. INTRODUCTION 



Supersymmetry (SUSY) is one of the most popular extensions of the Standard Model 
(SM) [l[ Q]. The main advantage with respect to non-SUSY models is the technical so- 
lution that it provides to the hierarchy problem [36]. In addition, other theoretical and 
phenomenological issues can be addressed in the context of supersymmetric models. For 
example, questions such as the radiative origin of the electroweak symmetry breaking [?-9] 



and the unification of the gauge couplings at high energies [10l4l4|. 



However, despite these appealing theoretical motivations, no experimental evidence of 
supersymmetry has been found so far at the Large Hadron Collider (LHC) jl^, Q. This 
might be telling us that we are adopting a wrong search strategy. In fact, most searches 
are based on the assumption of conserved R-parity [17|, thus including a lower cut on the 



amount of missing transverse energy. This should encourage the search for non-minimal 
supersymmetric scenarios with a departure from the usual signatures, such as those with 
R-parity violation [18 ]. 

In the Minimal Supersymmetric Standard Model (MSSM) there are renormalizable terms, 
perfectly allowed under all symmetries of the theory, that break either baryon or lepton 
number These new interactions, that involve SM particles and their superpartners, are 
highly constrained by the non-observation of baryon or lepton number violating processes 1 . 
In fact, if baryon and lepton number breaking interactions were present at the same time, 
they would induce proton decay. For these reasons, these dangerous terms are usually 
forbidden by an ad- hoc additional symmetry, such as R-parity. 

However, in order to stabilize the proton it is sufficient to forbid one of these terms, namely 
the baryon number violating one. In this case, the lepton number violating interactions 



an open issue 



would contribute to the generation of Majorana neutrino masses 18|, 
not addressed in the MSSM, while problems due to proton decay are evaded. This well 
motivated framework is one of the few scenarios where the mechanism behind neutrino mass 



generation can be directly tested in accelerator experiments, see for example 



29 



Furthermore, the presence of R-parity violating couplings opens up good perspectives for 



1 Many detailed studies regarding bounds on R-parity violating couplings have been written. These are 
obtained from baryon or lepton number violating processes as well as from several (unobserved) flavor 
violating processes. See for example 19-23 



25 27] and references therein. 
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a richer phenomenology. Although the smallness of the R-parity violating couplings implies 
that the standard decay chains are not affected, the decay of the LSP changes dramatically 
the final signatures at the LHC {37)]. In fact, the LSP decay reduces the amount of missing 
transverse energy in the final state. Instead, one expects events with large lepton and/or jet 



multiplicity. This has been recently used by different authors [35|, |38| to relax the stringent 
bounds on the squark and gluino masses, otherwise pushed to values clearly above the TeV. 

Regarding neutrino oscillation experiments, the intense activity over the last years has led 
to an increasing accuracy in the determination of the oscillation parameters. In particular, 
the #13 mixing angle has been finally measured, with a surprisingly large result. Indeed, 



Double-Chooz 



T2K 



4(j, MINOS |4l|, Daya-Bay |42j and RENO |43| have completely 



ruled out the possibility of a vanishing 6*13, consistently pointing towards a value in the 



sin 2 2#i3 ~ 0(0.1) ba 



appeared recently 



44 



lpark. Updated global fits to all available experimental data have also 



46]. 



The aforementioned measurement of #13 has also ruled out the popular tri-bimaximal neu- 
trino mixing pattern 47J and led to an explosion of papers where different flavor symmetries 
are used to accommodate a large #13. Concerning supersymmetric models with R-parity vio- 
lation, a flavor model for bilinear R-parity violation was introduced in 48(. The requirement 
of an exact A4 flavor symmetry in the superpotential couplings led to the known neutrino 
mixing pattern at the prize of a higher complexity of the scalar sector of the model. Here we 
present another example, based on a simple Z3 symmetry, of a flavored R-parity violating 
model. In addition to the bilinear term, the trilinear ones play a fundamental role in the 
generation of neutrino masses. The flavor symmetry strongly restricts the allowed terms, 
and only a few remain in the model. However, it still allows for the required interplay that 
leads to the observed neutrino masses and mixings. 

The paper is organized as follows. In section [Til we define the model and the basic 
assumptions followed along the paper. In section II III we describe the origin of the different 
contributions to the neutrino mass matrix and we analyze the resulting neutrino mixing 
pattern in section IIVI Finally, we briefly comment on collider phenomenology in section |V] 
and conclude in section [VT] 
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II. DEFINITION OF THE MODEL 



We consider the MSSM particle content with the following charge assignment under an 
additional Z 3 symmetry: 
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U 


Ei 




E3 


H d 


H u 


^3 
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id 


Lu 


1 


1 



TABLE I: Charge assignment of the model. We use the standard notation oj = e l2,r / 3 . 

In addition, all quark superfields are singlets under Z 3 . The most general superpotential 
in the leptonic sector allowed by the gauge and flavor symmetries is 

W = V/ /., /•.', //,/ + eLiH u + A, ; 7, 1., I. j !■'.!,■ + \'^ k LxQjD k , (I) 

where the only non-zero Xijk couplings allowed by Z 3 are A122, A133 and A231. The flavor 
symmetry also implies that only one e parameter is non-zero, e = e±, the other two being 
forbidden due to the \i and r charges under Z 3 . We do not specify the superpotential terms 
without lepton fields for brevity, since they do not play any role in the following discussion. 
We also impose the conservation of baryon number. This forbids the superpotential term 
X'-j k UiDjDk and estabilizes the proton. 

We allow for a soft breaking of the Z 3 flavor symmetry in the scalar potential, and thus 
general terms of the type 

Vsoft - rn 2 LiHd L*H d + B%H U1 (2) 

are considered. 

Finally, in the following discussion we will make two important assumptions about the 
parameters of the model. First, we will assume that the e parameter is a small dimensionful 
parameter and thus it plays a negligible role in neutrino mass generation. That will allow us 
to drop some contributions to the neutrino mass matrix, as required to obtain the correct 
neutrino mixing pattern. In the next section we will estimate how small e must be for our 
discussion to be unaltered. Second, we will consider universal soft SUSY breaking terms 
at the grand unification scale, incur- This implies that generation dependent couplings 
are exactly the same for all families and, in particular, one has rn 2 Lxlld = m 1 L2Ud = rr? LzHd 
and B\ = = B z e at m GUT . However, this equality is no longer valid at the SUSY scale 
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due to the effect of Renormalization Group Equations (RGE) running. In fact, in CMSSM- 
like scenarios, like the one considered here, 3rd family sfermions soft terms depart from 
universality at low energies due to their larger Yukawa couplings 2 . Therefore, we will simply 
make use of this departure from universality and consider independent 3rd family sfermion 
soft terms at the SUSY scale 3 . 



III. NEUTRINO MASS GENERATION 

The breaking of lepton number induces non-zero Majorana masses for the neutrinos 
which, in this model, are generated by both bilinear and trilinear R-parity violating cou- 
plings. As a result of that, the neutrino mixing pattern requires the correct interplay between 
those two sources. In this section we will describe the different contributions to the neutrino 
mass matrix and estimate their sizes. 

The interplay between different Ji p couplings is not a new idea. For example, a similar 
approach was followed in Ref . |3] , where the combination of bilinear and trilinear R-parity 
violation, together with a properly chosen flavor symmetry, also led to the required structure 
for the neutrino mixing matrix. However, the construction of the model and the resulting 
neutrino mass matrix are different in our work. The use of flavor symmetries to restrict the 



allowed tfl p couplings has also been discussed in [53l-l56|. 

Finally, before we proceed to describe the different contributions to the neutrino mass 
matrix we note that the charge assignment in table [J leads to a diagonal charged lepton 
mass matrix. 



A. Bilinear R-parity violation 

A precise determination of the 1-loop corrections to neutrino masses in bilinear R-parity 
violating supersymmetry requires the consideration of several types of diagrams. We will 



2 Another possibility that takes non-universality in the soft terms as starting point is the so-called natural 
supersymmetry, where the 1st and 2nd generation sfermions are much heavier than the 3rd one, see for 



49 



50 



example 

3 Approximate formulas for the 3rd family soft terms are known in the literature, see 51] . However, their 



complexity (and dependence on unknown parameters such as tan 0) does not allow for simple estimates. 
Therefore, we will treat the departure from universality as a free phenomenological parameter. 
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FIG. 1: 1-loop neutrino masses from bilinear R-parity violation. B e contribution. 
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FIG. 2: 1-loop neutrino masses from bilinear R-parity violation. m 2 LH contribution. 



follow closely the results presented in the pioneer references on this topic |57H60|. 

The Z3 flavor symmetry implies that only one bilinear e parameter is allowed in the 
superpotential, eL\H u . This would in principle lead to one non-zero entry in the tree- level 
neutrino mass matrix, m^ 1 oc e 2 . However, we will assume that the smallness of the e 
parameter makes this contribution negligible. 

Therefore, the dominant BRpV contributions to neutrino masses come from 1-loop di- 



agrams without the e coupling (see Ref. 6l| for a similar scenario). Examples of such 



possibilities are presented in figures [T] and [2j We note however that they require relatively 
large B\ , vn 2 L _ Hd (see below for an estimate), otherwise these 1-loop contributions would 
be too small to account for neutrino masses. This would generate large sneutrino vacuum 
expectation values (VEVs), v l L , and lead to too large tree-level contributions to neutrino 
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masses. Therefore we are forced to impose a fine-tuning of the parameters so that the rel- 
ative contributions to v l L coming from B\ and m? L . H cancel each other. This fine-tuning 
relation can be easily obtained from the potential in equation (J2J), 



ml iHd v d + B\v u ~ 0, 



(3) 



which implies that the approximate relation m 2 LiHd ~ tan f3 B\ must hold. 

We proceed now to estimate the size of the different diagrams. If the scale of lepton 



number violation in the soft terms is denoted as B\ , rn 2 LHd 



m 



\ , one obtains 4 



m 



BRpV, 



1-loop) 



m, 



16tt 2 m| USY ' 



(4) 



For msusY — 100 GeV one finds that m# ~ 1 GeV leads to m BRpV ~ 0.1 eV, in the correct 
range. On the other hand, the tree-level contribution to the neutrino mass matrix can be 
roughly estimated to be m BRpV (tree) ~ e 2 /msusY, and thus for msusY = 100 GeV one needs 
e > 100 KeV to be competitive with the 1-loop contributions. We will therefore assume in 
the following that e is below that scale. 

It is important to notice that these contributions require a non-vanishing mass splitting; 
between the real and imaginary components of the sneutrinos, Am?, as discussed in [57J and 
previously deduced from general considerations in 621 ]. 

Concerning the flavor structure of these 1-loop contributions, note that by assumption 
the soft SUSY breaking potential is flavor universal for the 1st and 2nd generations, with 
a (little) departure in the 3rd generation soft terms. Therefore, we obtain the sum of two 
terms, a democratic structure plus a deviation from universality 



ni 



BRpV 



(l i x\ 

1 1 1 
1 1 1 



d 



( 



\ 



1 \ 
1 
112; 



(5) 



Here a is given by the estimate in equation (j3J). The dimensionful coefficient d follows from a 
similar expression, where the corresponding soft terms m 2 ^ have been properly replaced by 
^ +5m 2 . The mass squared parameter 5m 2 is the result of the deviation from universality, 



m 



4 In fact, this definition of rnL^ actually corresponds to m 2 LiHd . That is why equation (TJ]) does not include 
a tan 2 ft enhancement, see 0, [|3(]], since that factor is absorbed in m 2 ^, . 
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expected to be sizable for the 3rd generation. Terms of order {8m 2 ) 2 have been neglected in 
equation (jSJ). 

Before we move on to the discussion of the contributions from trilinear R-parity violating 
couplings, we would like to emphasize that the deviation from universality provided by the 
second term in Eq. (J5]) is phenomenologically required. As we will see below, this is the 
only term that breaks \x — r invariance in the neutrino sector. Its absence would imply a 



vanishing reactor mixing angle and maximal atmospheric mixing, see for instance [63|]. This 
would be in contradiction with recent data, that clearly indicates #13 7^ 0. Therefore, one can 
interpret the departure from universality in the soft terms, quite general in supersymmetric 
models, as the origin of the non-zero reactor angle in our model 5 . 



B. Trilinear R-parity violation 

As previously explained, the flavor symmetry implies that the only non-zero Xijk couplings 
are A122, A133 and A231. Therefore, the usual 1-loop diagrams in trilinear R-parity violation, 
see figure [3l can only fill the 11 and 23 — 32 entries of the neutrino mass matrix, leading to 

( h 0^ 



m TRpV 



0c 
yO c j 



(6) 



Similarly, the A' couplings can also contribute to neutrino masses, but only to the 11 
element of the matrix, here denoted by b. We can now estimate the size of this family of 
1-loop contributions. For example, in the case of the c element in the previous matrix one 
has 

K RPV ) 23 = c ~lh A^ tan/3 Ht^L. (7) 

The element b has a similar generic expression, although more couplings are involved (and 
for the A' diagrams quark masses appear instead). Equation © shows that 0(0.1 eV) 



5 As already explained, the fi — t invariance can also be broken by hand by introducing a small difference 
between the soft SUSY breaking terms B"j and flf (or, analogously, Tn\^ n an< ^ m l 3 j? d )- Similarly, one 
can relax the conditions given by the Z% symmetry by allowing the existence of T\ couplings whose 
corresponding superpotential terms are forbidden. For example, T\ 233 would induce a high order breaking 
of the fJ> — t invariance of the neutrino mass matrix. However, we consider the breaking by the Yukawa 
couplings, transferred to the soft SUSY breaking terms by RGE running, the simplest solution. 
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FIG. 3: 1-loop neutrino masses from trilinear R-parity violation. 

contributions can be obtained for /i ~ msusY = 100 GeV and tan/3 = 10 if A231 ~ 0.05. 
This value is close to the experimental limit, but it can be easily lowered by using \i > m SU sY 
or larger tan/3. See table HTT1 and refs. 2l|, |23( for more details. 
Summing up the two contributions one obtains the texture 



m BR P V + m TRpV 



a + b a a + d 

a a a + c + d 

a + da + c + d a + 2d 



\ 



(8) 



In the following section we will show how the mass matrix in Eq. (JS]) can accommodate 
the observed pattern of neutrino mixing. However, we already observe that the coefficient 
d, which can be traced back to the departure from universality in the soft SUSY breaking 
terms, is the only source of fi — r invariance breaking. In fact, in the limit d = 0, one 
recovers an exactly fi — r symmetric neutrino mass matrix that implies #13 and maximal 
atmospheric angle 6 . In conclusion, d 7^ is required. This will induce a non-zero reactor 
angle and a deviation from maximal atmospheric mixing. Since the same parameter is at 
the root of both deviations, we expect a tight correlation between #13 and #23- 



We remind the reader that we are working in a basis where the charged lepton sector is diagonal. 
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IV. NON-ZERO 0i3 AND DEVIATIONS FROM MAXIMAL ATMOSPHERIC 
MIXING 

The mass matrix in Eq. ([S]) has only four parameters. Therefore, since it must be able to 
accommodate the pattern of neutrino mixing measured in oscillation experiments, some con- 
nections between the parameters must appear. Although we cannot give definite predictions 
for the oscillation parameters, clear correlations between them exist. Similarly, the absolute 
scale of neutrino masses, determinant for neutrinoless double beta decay experiments, is 
linked to the other parameters. This allows to test the model with neutrino phenomenology. 

In our analysis we assumed real parameters 7 , so the number of free parameters in the 
neutrino sector is 4, namely, a, b, c, and d in Eq. ([H]), to explain 6 observables, namely two 
squared mass differences Am 2 im and Am^, the neutrinoless double beta decay effective 
mass m ee , and the three mixing angles. Therefore we have two predictions: the neutrinoless 
double beta decay effective mass and a correlation between the reactor mixing angle and 
the atmospheric mixing angle. 

We have performed a detailed scan over the parameter space fixing Am^ m , Am^, sin 2 9u 



and sin 2 # 2 3 within the 3a range given in [44|. Large deviations from the measured value of 
#13 have been allowed in order to clearly see potential correlations between this parameter 
and the other mixing angles. Correct fits 8 have been found for both, normal and inverted 
neutrino spectra, without a particular preference. However, the results clearly favor a quasi- 
degenerate spectrum with large neutrino masses. 

The mass matrix in Eq. ([8]) would be /i — r symmetric in the absence of the d parameter, 
thus leading to sin 2 6*23 = 1/2 and sin 2 #13 = with a free 6*12. Interestingly, deviations 
from this scheme typically affect simultaneously #13, that departs from zero, and #23, that 
departs from maximality. This generic prediction, present in many flavor models (see for 
instance [24], is especially relevant in our model, since it is caused by a single parameter, 
d, which leads to a tight correlation. Our numerical results confirm this expectation, see 
Fig. HI where sin 2 #23 is shown as a function of sin 2 #13 for normal (left-side) and inverted 



7 This assumption will be lifted below when the general CP violating case is considered. 

8 In the following, we will use the expressions valid fit or correct fit to denote a set of values for the 
parameters a, b, c and d that leads to Am 2 tm , Am^,, sin 2 # 12 and sin 2 2 3 within the 3a range given in 



44] 
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FIG. 4: sin 2 ^23 as a function of sin 2 #13 for normal (left-side) and inverted (ri ght- side) neutrino 



4J. The dashed 



spectra. The numerical scan is based on the results of the global fit of Ref. 
lines represent the best-fit values for sin 2 #13 and sin 2 #23 1 whereas the shaded regions correspond 
to 3a (in case of sin 2 #13) and la (in case of sin 2 #23) allowed regions. Different colors correspond 
to different CP branches: r\\ (blue), 772 (red) and 773 (purple). 

(right-side) neutrino spectra. The atmospheric mixing angle is restricted to be within the 
3cr allowed region found in [441 ] . and thus sin 2 #23 £ [0.39,0.64], whereas the vertical shaded 



band represents the 3<r allowed region for sin 2 #13, as given by the same analysis. 

Figure H] deserves some comments. In addition to the aforementioned correlation, one 
clearly observes distinct regions in the #13 — #23 plane. These re gion s, or branches, correspond 



to different cases of intrinsic CP charges of the neutrinos, 77 (64j. There are four possible 
cases: r]i = (-, +, +), rj 2 = (+, — , +), ^3 = (+, +, — ) and 774 = (+, +, +). Other choices can 
be reduced to one of these by means of an unphysical global sign. In Fig. HI different colors 
are given for these CP cases: blue for 771, red for 772 and purple for 7/3. No valid fits were 
found for the 774 case. Similarly, for inverted hierarchy we could only find correct points for 
cases 771 and 772. 

For normal hierarchy, the 7/4 case leads to very large deviations in #23 that clearly departs 
from maximality even for small values for 6* 13 . In fact, from our results in Fig. H]one can 
conclude that the 771 branch is ruled out by oscillation data since it cannot accommodate 
6*i3 and #23 simultaneously. In contrast, cases r] 2 and 773 can reproduce the measured mixing 
angles, although with different predictions. On the one hand, 772 leads to large deviations 
in #23, which must lie 3a away from the best-fit value in order to be compatible with 
the measured value for # 13 . On the other hand, 773 leads to small deviations, at most 
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Analysis 



Best-fit ±la 3cj range 



Forero et al. [44 



0.49 



0.53 



+0.08 
-0.05 

-0.05 
-0.07 



0.39 - 0.64 



Fogli et al. [45] 



- 386 ±o:o2i °- 331 - °- 637 



0.392 



-0.039 
-0.022 



Gonzalez-Garcfa et al. [4g] 0.42+°;°^ 



0.335 - 0.663 
0.39 - 0.66 



TABL 
et al. 

to inverse hierarchy. 



II: Comparison between the results of Refs. 
1 

44] and Fogli et al. 



44 



46 ] regarding sin 2 623 ■ In case of Forero 



451 ] the upper row corresponds to normal hierarchy and the lower row 
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FIG. 5: sin 2 #23 as a function of sin 2 #13 for normal (left-side) and inverted (ri ght 



side) neutrino 



spectra. The numerical scan is based on the results of the global fit of Ref. [45J. The dashed 
lines represent the best-fit values for sin 2 #13 and sin 2 #23 > whereas the shaded regions correspond 
to 3(7 (in case of sin 2 #13) and la (in case of sin 2 #23) allowed regions. Different colors correspond 
to different CP branches: r\\ (blue), 772 (red) and 773 (purple). 

A sin 2 #23 ~ 0.2. For inverted hierarchy the results are slightly different. First, it is now 772 
the one that is ruled out due to the too large deviations induced in Q23 ■ And second, for rji 
one finds very small departures from maximal atmospheric mixing, at most A sin 2 623 ~ 0.03. 

So far we have based our discussion on the results of the global fit of Ref. 44j . In their 
analysis, the authors of this work found no significant deviation from maximal atmospheric 
mixing, with best-fit values sin 2 623 = 0.49 (for normal hierarchy) and sin 2 623 = 0.53 (for 



inverse hierarchy). However, other recent global fits 45|, |46| find a slight preference for 
#23 < tt/4. For example, in Ref. 45[ the results of the analysis show a < 2a deviation from 
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maximal mixing in case of normal hierarchy, with best-fit value sin 2 #23 = 0.386, and a < 3<r 
deviation in case of inverse hierarchy, with best-fit value sin 2 #23 = 0.392. Similarly, in the 



preliminary results presented in 46| one can also find a preference for #23 < 7r /4, although 
the statistical significance is lower. A brief compilation of these results is given in Table [TT1 
Due to the impact that such deviation from maximal atmospheric mixing would have on 
our model, because of the correlation with the reactor angle, we have repeated our numerical 



scans using the input parameters given by Refs. [45|, |46|. The results for the case of Ref. (45| 
are shown in Fig. [5j Qualitatively the results do not differ very much from those presented 
in Fig. HJ There are, however, some differences. For example, there is a little preference for 
the 773 case and a neutrino spectrum with normal hierarchy. We do not show the analogous 
figures for Ref. {46], since they lead to very similar conclusions. 

We have also investigated the predictions that our model can make regarding neutrinoless 
double beta decay (0u2/3). The results are shown in Fig. [6j where the effective neutrinoless 
double beta decay parameter m ee is given as a function of the lightest neutrino mass. We 
present the results for normal hierarchy on the left-side and the results for inverse hierarchy 
on the right-side. The shaded bands (in grey and yellow) correspond to the well-known 
flavor-generic predictions, much more spread than the predictions given by our model. In 
fact, one can easily see that our model is restricted to lie on a very small portion of the 
plane. This is due to the fact that only solutions with quasi-degenerate spectra were found 
in our scan. 

Again, different CP intrinsic charges for the neutrinos lead to different regions in Fig. |6j 
The color code is as in Fig. HI Future experimental sensitivities are also included in this 
figure, as well as the region excluded by the Heidelberg-Moscow collaboration [65|. Since our 
model favors quasi-degenerate spectra, most of the points have large values for the lightest 
neutrino mass and the effective rn P , P Ou20 parameter. This implies that they are at the reach 



of experiments such as KATRIN [60], GERDA-II [67j and EXO-200 



68| 9 . Only 7/3 in the case 



of normal hierarchy predicts low values for m ee , out of reach for GERDA-II and EXO-200, 



9 Recently, the EXO-200 collaboration made public new results 69(, not included in Fig. [BJ which allow to 
set a lower limit on the 0^2/3 half-time in 136 Xe of about T 1 / 2 > 1.6 x 10 25 yr (90 % CL). This corresponds 
to m ee < 140 — 380 meV, the broad range being caused by the theoretical uncertainty in the nuclear 
matrix element calculation. Similarly, the KamLAND-Zen collaboration also reported new limits [toI | . 
although slightly less stringent. 
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FIG. 6: Effective neutrinoless double beta decay parameter function of the lightest neutrino 

mass for normal (left-side) and inverted (right-side) neutrino spectra. The grey and yellow shaded 
regions correspond to the flavor-generic normal and inverse (gray) hierarchy neutrino spectra, 
respectively. For a discussion see text. 

but still with sufficiently large neutrino masses so that KATRIN is able to measure them. 
In conclusion, both Oz/2/3 and absolute neutrino mass experiments offer good perspectives 
to probe our model. 

Finally, let us briefly discuss the general CP violating case, where all four parameters in 
the neutrino mass matrix are allowed to be complex. We have generalized our numerical 
scan to include all CP phases, the Dirac, 5, and Majorana ones, a and /3. The main result is 
that the regions between branches in Figs. H] and [5] are filled with points where at least one 
of the (physical) CP phases do not vanish. Therefore, the concluions drawn in the previous 
paragraphs should be regarded as valid only in the CP conserving case. For example, our 
model would not be ruled out if the #i 3 and 9 2 3 angles are precisely measured to lie outside 
the 771,2,3 branches. That would point towards a CP violating scenario, something to be 
probed by experiments such as T2K 71 1 and NOz/A [^j], see \ 
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Coupling 


Bound 


The bound comes from 




2.7- 10" 2 


Neutrino masses [211 


A133 


1.6 • 10~ 3 


Neutrino masses [211 


A231 


0.05 


Leptonic r decay [23] 


Xfk (j + k) 


0.02 - 0.47 


Several processes [21] 


* 3 k U = k) 


3.3 • 10~ 4 - 0.02 


fifi0i> and m u [21] 



TABLE III: Experimental bounds for trilinear R-parity violating couplings. These limits were 
obtained by setting all SUSY masses to 100 GeV. For the A' couplings the extreme cases are 
shown, with the most and least constrained couplings for each case. For more details see references 
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23[]. 



V. SOME COMMENTS ON COLLIDER PHENOMENOLOGY 

Since R-parity is broken the LSP is no longer stable and decays 10 . The decays can go via 
bilinear (e, B\ and m 2 L _ Hd ) or trilinear (Ai 22 , A133, A231 or A^) fi p couplings. In the bilinear 
case, only e can mediate the decay at tree-level, and thus its small size reduces the relevance 
of this possibility. We are thus left with decays mediated by trilinear couplings and we 
conclude that gauge-mediated LSP decays are suppressed. 

By studying the size of the different A and A' couplings one could go further and find 
the dominant LSP decay signatures that particular scenarios would predict. For example, 
a simple estimate based on the 1-loop generated masses that come from A231 shows that 
this coupling must be of order 0.05. On the other hand, the other two A parameters (Ai2 2 



and A133) are strongly constrained by flavor physics 21], and thus they are less relevant for 



the LSP decay. However, the relative importance of A231 and the X' jk couplings cannot be 



determined. Tab 
references 
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e lllll shows the relevant experimental bounds for our model, as obtained in 
23j. However, we cannot know a priori the relative importance of the different 
couplings and general predictions concerning the LSP decay cannot be made (apart from 
the aforementioned trilinear dominance). 

Nevertheless, we emphasize that independent tests of the model and its underlying fla- 



10 This implies that the usual neutralino LSP is lost as dark matter candidate. However, it is well-known 
that a relatively light gravitino provides a valid alternative, see for example the recent references 74-76$. 
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vor structure should be performed. In addition to the aforementioned correlations among 
neutrino mixing angles, unfortunately present in many flavor models, collider tests are fun- 
damental in order to disentangle the dynamics behind the observed flavor pattern in neutrino 
mixing. Reference addresses this issue by studying how one may discriminate between 
different $ p flavor operators leadin g to the decays of a neutralino LSP at the LHC. Simi- 
lar works exist in the case of b-fc [23-133, y, l36|, where one can find a strong correlation 



between LSP decays and neutrino mixing angles. Finally, it is also possible to distinguish 
between bilinear and trilinear violation of R-parity if the LSP is a slepton. In that case, its 
decays properties can be used to determine the relative importance of the tfL p couplings, as 



shown in 



33]. 



VI. CONCLUSIONS 



We studied a simple extension of the MSSM assuming an abelian Z% flavor symmetry. 
Neutrino masses arise at 1-loop from bilinear and trilinear R-parity breaking operators. 
If the soft SUSY breaking terms were flavor blind, one would obtain a // — r invariant 
neutrino mass structure, thus implying a vanishing reactor angle. However, this is cured 
by renormalization group running effects, naturally present in supersymmetric theories. We 
have shown how this model can accommodate a neutrino mixing pattern in accordance to 
data. In particular, we have studied the correlations between mixing angles that allow us 
to link different deviations from tri-bimaximal mixing. The model also predicts large values 
for the lightest neutrino mass and the effective m ee Oz/2/3 parameter, both at the reach of 
current and future experiments. Finally, since the LSP decays one in principle has a rich 
phenomenology at the LHC. Due to the structure of the model, we predict a clear dominance 
of the trilinear couplings in the LSP decay. 
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